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We investigate the bulk magnetic, electron paramagnetic resonance, and magneto-optical prop- 
erties of {Ni4Moi2j-, a magnetic molecule with antiferromagnetically coupled tetrahedral Ni 11 in a 
diamagnetic molybdenum matrix. The low-temperature magnetization exhibits steps at irregular 
field intervals, a result that cannot be explained using a Heisenberg model even if it is augmented 
by magnetic anisotropy and biquadratic terms. Allowing the exchange and anisotropy parameters 
to depend on the magnetic field provides the best fit to our data, suggesting that the molecular 
structure (and thus the interactions between spins) may be changing with applied magnetic field. 
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I. INTRODUCTION 



Despite the enormous progress that has been made 
in understanding magnetic molecules over the past 
decade j b 2 ! 3 ! 4 ! 5 ! 6 ! 7 ' 8 ' 9 it is still a challenge to deduce the 
underlying microscopic spin Hamiltonian. Mni2-acetate 
is a good example of this problemi^iMiiiSiH Although 
known for almost twenty years, only extensive investi- 
gation elucidated the model parameters Small mag- 
netic molecules, with their simpler chemical and mag- 
netic structure, can therefore provide an important op- 
portunity to understand the dependence of magnetic ob- 
servables on model parameters ^LiLi2i 

In this work we report our joint experimental 
and theoretical efforts to understand the behavior 
of [Moy 2 03o(At2-OH) 10 H 2 {Ni n (H 2 0)3}4], henceforth ab- 
breviated as {N14M012}, a magnetic molecule which is 
comprised of Ni 11 centers positioned at the nucleophilic 
sites of an e-Keggin cluster forming an almost ideal 
tetrahedron. 19 In contrast to several other nickel com- 
pounds which exhibit ferromagnetioi^^SiSiiSSiS 3 " or mixed 
coupling^ the Ni centers of this molecule are anti- 
ferromagnetically coupled, as is also the case for cer- 
tain Ni-2 x 2-grid molecules^ Because the structure of 
{N14M012} is almost perfectly tetrahedral one might an- 
ticipate that the magnetic energy levels are reasonably 



well described by an isotropic Heisenberg model with a 
single antiferromagnetic exchange parameter^ Such a 
Hamiltonian can be expressed in terms of the total spin: 



H = -2J 
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where s(u) is a single-spin operator at site u and S is 
the total spin operator. The spin quantum number of 
each Ni 11 ion in an octahedral ligand field is s = 1. 
For antiferromagnetic coupling (J < 0) the resulting 
low-tempcrature magnetization curve M. (B) that follows 
from displays four steps before reaching saturation. 
These steps occur at the magnetic fields 
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for S — 0,1,2,3, where the lowest Zeeman-split lev- 
els of adjacent multiplets cross. In particular, from 
(J2J it follows that the level crossing fields are uni- 
formly spaced. 25 26 In stark contrast to the expecta- 
tion of uniformly spaced crossing fields, the experimental 
magnetization A4(B) curve of {M4M012} features non- 
cquidistant steps at 4.5, 8.9, 20.1, and 32 T. Even as- 
suming an anisotropic Hamiltonian^ 7 . with two exchange 
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couplings and biquadratic terms as done for other Ni 11 - 
compoundsji£ii& we were unable to account for the spe- 
cific sequence of steps in the low-temperature magnetiza- 
tion of {M4M012}. In order to provide a comprehensive 
picture of the unusual high-field magnetic behavior of 
{M4M012}, we have been led to invoke field-dependent 
exchange and anisotropy parameters. We argue that this 
dependence emanates from changes in molecular struc- 
ture (and thus the overlap of those atomic orbitals that 
determine the exchange interactions as well as the coor- 
dination geometries which affect the electronic single-ion 
properties) with applied magnetic field. The magneto- 
optical response of {M4M012} supports a small change 
in the Ni 11 coordination environment and the associated 
electronic single-ion properties. 



II. CHEMICAL AND CRYSTAL STRUCTURE 




FIG. 1: (color online) Ball-and-stick representation of the 
{Ni4Moi2 } molecule (Ni: numbered big spheres; Mo: medium 
size spheres; O: small spheres; H: not shown) emphasizing a 
slightly stretched Ni 4 : pyramid with a near-equilateral trian- 
gle base (Ni ions labeled 2, 3, and 4) and an elevated apex 
(Ni ion labeled 1). Ni-Ni distances: d 12 = 6.700(5) A, d 13 = 
di4 = 6.689(1) A, d 23 = d 24 = 6.616(1) A, d 34 = 6.604(1) Ap±2 



The neutral {N14M012} cluster is isolated in the form of 
crystals of [Moy 2 3 o(M2-OH) 10 H 2 {Ni n (H 2 0)3}4]-14 H 2 
and is based on the diamagnetic, highly-charged e-Keggin 
anion [Mo^Osg^s-OH)^ 18 ", built up from four edge- 
sharing {M03} groups (each consisting of three edge- 
sharing MoC>6 octahedra). Within the e-Keggin frame- 
work, the Mo positions form six Mo^ groups with short 
Mo-Mo single bonds. The e-Keggin structure is formally 
derived from the common a-Keggin isomer by rotating 
all four M03 groups by 60°, preserving the T^ sym- 
metry. In {N14M012}, four [Ni n (H 2 0)3] 2+ groups are 
coordinated each to three (unprotonated) /x 2 -oxo cen- 
ters that interlink the Mo positions of the Mo^ groups 
(Fig. ID . This results in a octahedral 3 Ni n (H 2 0)3 coor- 
dination environment with all-trans-positioned oxo and 



water ligands and nearly identical Ni-0 distances (Ni- 
(M2-O): 2.05 A, Ni-OH 2 : 2.06 A). This capping of the 
M03 groups of the e-Keggin fragment by four Ni 11 po- 
sitions produces a near-regular N14 tetrahedron. Con- 
trary to many other tetrahedral M4 structures in which 
Ni pairs are connected by mononuclear bridging centers, 
in {Ni4Moi 2 } each Ni pair is interconnected via one - 
0(-Mo-) 2 0- bridging motif serving as a superexchange 
pathway (see Fig. ^ where one Ni-O-Mo-O-Ni pathway 
is highlighted by dark bonds). The geometry of each of 
these pathways is therefore characterized by four bond 
lengths, three bond angles, and two dihedral angles, as 
opposed to Ni4-type structures comprising mononuclear 
linker groups (two bond lengths, one bond angle). Im- 
portantly however, the molecular geometry of {Ni4Moi 2 } 
slightly deviates from the T^-symmetric ideal, resulting in 
a slightly stretched M4 pyramid with elongated Ni-Ni dis- 
tances between the Ni positions of a basal M3 plane and 
one apex. While the crystallographic symmetry opera- 
tions result in the molecular point group C s , the actual 
geometry is virtually of C% v symmetry, and the geomet- 
ric parameters for all Ni-0(-Mo-) 2 0-Ni pathways corre- 
spondingly fall into two sets. Within these two sets, the 
individual bond lengths and angles display minimal devi- 
ations (typically < 0.8 %) from the respective averages: 
Intra-basal Ni-Ni contacts are characterized by (Ni-O) = 
2.05A, (Mo-O) = 1.95 A, (Ni-O-Mo) = 135.4°, and (Mo- 
O-Mo) = 89.6°. For Ni-Ni contacts between the apex 
to the base positions average values of (Ni-O) = 2.03 A, 
(Mo-O) = 1.95 A, (Ni-O-Mo) = 137.9°, and(jMo-O-Mo) 
= 94.6° are found. As the geometric parameters do not 
vary significantly within the intra-basal and within the 
apex-basal Ni-Ni contacts, we do not take into account 
slight deviations from the idealized C^ v symmetry, and 
use only exchange constants J and J' in this paper (see 
Fig. . 

,f. 1 
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FIG. 2: Simplified structure of Ni4: the superexchange inter- 
actions J' and J are represented by dashed and solid lines, 
respectively. 



It should be added that magnetic exchange through 
polyoxomolybdates frameworks, especially if more than 
one possible pathway exists for each contact, has been 
found to be fairly insensitive to the separation distance of 
the pair of spin centers. Also, for similar systems based 
on mononuclear linker groups the bond angles have a 
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stronger effect on the exchange energies than the con- 
tact distance. 16 i 18 i 20 i 21 i 22 i 28 i 29 i 30 Due to the presence of 
crystal water molecules in the solid-state structure of 
{Ni4Moi2}-14 H2O which space the cluster entities apart, 
the closest intermolecular Ni- • • Ni distance in the solid 
state exceeds 7.15 A, rendering inter-molecular (dipolc- 
dipole) magnetic exchange insignificant. 

III. EXPERIMENTAL RESULTS 
A. Magnetic properties of {NuMo^} 
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FIG. 3: (color online) Low-field susceptibility: Experimen- 
tal data are given by black squares whereas the solid curve 
shows the result assuming a simple Heisenberg Hamiltonian 
Q with J/fcs = —3.4 K and an isotropic spectroscopic split- 
ting factor g = 2.25. The inset shows a close-up view of the 
low-temperature region for different magnetic fields. 

Figure displays the magnetic susceptibility M(B) 
of {M4M012} as a function of temperature. These data 
were collected on a SQUID magnetometer (Quantum De- 
sign MPMS-5) at various magnetic fields in a tempera- 
ture range of 2 — 290 K. Using a simple Heisenberg Hamil- 
tonian Q one obtains J/ks = —3.4 K and an isotropic 
g-factor of g — 2.25, see also Ref. 19J The model fit given 
by the solid curve in Fig. [31 is acceptable, although the 
low-temperature behavior is not well reproduced, as can 
be seen in the inset of Fig. |3 where the magnetic suscep- 
tibility for various field values and low temperatures is 
displayed. We suggest that the deviation is not only due 
to the presence of impurities but also due to anisotropic 
as well as biquadratic terms in the Hamiltonian which 
are known to be needed to model the low-temperature 
behavior of Ni II -compoundsiiSii 8 . 

High-field magnetization measurements are a valuable 
tool to extract information on the spin Hamiltonian 
which is not accessible from magnetization measurements 
on commercial SQUID magnetometers. The high-field 
magnetization for a powder sample of {M4M012} has 
been independently measured in pulsed magnetic fields at 
the facility of the National High Magnetic Field Labora- 
tory (NHMFL) at Los Alamos as well as at the Okayama 
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FIG. 4: (color online) Magnetization: Experimental data are 
given by squares (dark symbols - NHMFL, T = 0.44 K; light 
symbols - OHMFL, T = 0.40 K). The theoretical magne- 
tization assuming a simple Heisenberg Hamiltonian Q with 
J/ks = —3.4 K and an isotropic spectroscopic splitting factor 
g — 2.25 is given by a solid curve for T = 0.44 K. The inset 
shows the experimental differential magnetization AM /AB as 
data points as well as the theoretical AA4/AB (solid curve) 
using the same parameters as above. 



High Magnetic Field Laboratory (OHMFL) by using a 
standard inductive method (maximum at NHMFL B = 
60 T, whereas the maximum at OHMFL is B = 40 T, 
dB/dt = 10000 . . . 15000 T/s). The results of these two 
measurements are in very close agreement. No hystere- 
sis is found between up and down sweep runs, indicating 
thermal equilibrium behavior of M. vs. B. 

Figure 0| shows the magnetization as a function of ap- 
plied external magnetic field at T = 0.44 K. Strikingly, 
four non-equidistant steps are observed in the magne- 
tization. These steps are found near 4.5, 8.9, 20.1, and 
32 T. Saturation of the magnetization is not observed un- 
til 60 T. For comparison, we show the expected response 
of a simple Heisenberg Hamiltonian (Eq. ideal tetra- 
hedron) for the model parameters extracted above. Note 
that this model predicts equidistant steps in the magneti- 
zation at 4.5, 9.0, 13.5, and 18.0 T with saturation above 
the fourth step. The drastic deviation between the pre- 
dictions of Eqn. and the experimental results cannot 
be the result of heating via the magnetocaloric effect. 
Although observed in other compounds, such an effect 
would only smear out the steps but not shift the step po- 
sitions. In addition sample heating or cooling in a varying 
field is often accompanied by hysteresis, which is absent 
in this measurement;— The observed step positions thus 
constitute a challenge not only to the simple Heisenberg 
model given by Q), but also to more elaborate mod- 
els © incorporating anisotropy terms and biquadratic 
exchange. We suggest that these models should be ex- 
tended to include field-dependent parameters in order to 
account for our experimental results. This will be dis- 
cussed in Sec. 
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B. EPR response of {Ni 4 Moi 2 } 
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FIG. 5: (color online) The figure shows the dependence of the 
observed EPR resonance frequencies PI (circles), P2 (trian- 
gles), and P3 (squares) on the magnetic field. The lines pro- 
vide linear fits to the data. The data were taken at T = 4.2 K. 

Figure [S] displays the results of our Electron Paramag- 
netic Resonance (EPR) measurements. The transmission 
of a powder sample was determined as a function of ap- 
plied magnetic field in a frequency range from 95 GHz 
to 381.5 GHz. The figure shows the dependence of the 
observed EPR resonance frequencies PI (circles), P2 (tri- 
angles), and P3 (squares) on magnetic field. One imme- 
diately notices that two different slopes can be assigned 
to the data, one corresponding to AM = 1 and another 
one corresponding to forbidden transitions with AM = 2. 
These dependencies can be qualitatively - and to some 
extent quantitatively - explained by looking at the Zee- 
man level scheme of the simple Heisenberg model as 
it is schematically depicted in Fig. HJ1 




FIG. 6: (color online) The figure shows schematically the Zee- 
man level splittings in a pure Heisenberg model together with 
the assignments of allowed and forbidden transitions. 

The strongest transition which we observe in the spec- 
tra is the allowed transition P3 with AM = 1. Since 
the applied temperature is of the order of the coupling, 
this transition is actually a sum of several transitions. At 
low-field values it is dominated by the transition between 



(S = 1, M = -1) and (S=1,M = 0), whereas at higher 
fields the dominant contribution stems from the transi- 
tion between (5 = 2,M = -2) and (5 = 2, M = -1). 
The dependence of the resonance frequency of P3 on the 
applied field suggests that the zero-field splitting in the 
triplet is small. 

PI is a low-field transition which connects (S = 1, M = 

— 1) and (S = 1, M = +1) and thus should be forbidden. 
In the spectra its strength is much weaker than that of 
P3. We believe that this transition appears due to mixing 
of S z eigenstates that would arise from anisotropic con- 
tributions to the Hamiltonian. The line which is plotted 
through the data points suggests that the zero-field split- 
ting between (S = 1,M = -1) and (S = 1,M = +1) is 
probably small, although, the experimental data points 

- which extend only down to 1.3 T - would also allow 
a somewhat bigger splitting, especially since the lowest- 
lying PI data points appear to deviate from a straight 
line. 

P2 is another rather weak forbidden transition which 
shares the slope with PI. We believe that this transi- 
tion connects (S = 2, M = -2) and (S = 1,M = 0). 
This transition is not observable below about 7 T due 
to the fact that this transition occurs only when the 
mixing of S z eigenstates is sufficiently strong which is 
the case around the level crossing at 9 T. The depen- 
dence of this transition on temperature, i.e. on the ther- 
mal occupation of the level with (S — 2,M = —2) is 
small. The extrapolation of the field dependence al- 
lows one to deduce an approximate isotropic Heisen- 
berg coupling from the zero-field energy separation of 
the triplet and the pentuplet. From Eq. one deduces 
that E(S = 2) - E(S = 1) = -4J, thus one obtains 
J w —3.4 K, which is in very good agreement with the 
exchange constants deduced from our susceptibility mea- 
surements. The spectroscopic splitting factor can be de- 
termined from the slopes of P3 and PI in Fig. to be 
g = 2.23 ±0.03, which is also in good agreement with the 
value deduced from the high-temperature behavior of the 
susceptibility. It is noteworthy that the high field data of 
transition P3 have a smaller slope which, if fitted alone 
(dashed-dotted line in Fig.EJ, suggests g w 2.11 ± 0.03. 
However, we want to point out that these considerations 
are done on the basis of the simple Heisenberg model Q), 
that does not take into account that a realistic Hamilto- 
nian has to contain anisotropic terms. A more detailed 
explanation is given at the end of sec. 

Summarizing this part, we find that the zero-field split- 
ting is small. We explain this observation by the fact that 
the local principal axes of the Ni-ions point in different di- 
rections (radially outwards), and thus the average global 
anisotropy is small. Our use of a powder sample led to 
additional averaging. 
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FIG. 7: (color online) Optical conductivity of a {N14M012} 
pellet at 4.2 K (solid curve) and 300 K (dashed curve), cal- 
culated from reflectance measurements (inset) by Kramers- 
Kronig analysis. The energy range of our magneto-optical 
work is indicated by the arrow. 



C. Optical properties of {N14M012} 

Figure [7| displays the optical conductivity of 
{NI4M012}. These experiments were carried out on a 
pressed powder sample using a Lambda-900 grating spec- 
trometer equipped with reflectance stage and cryostat^i 
{M4M012} is a semiconductor with an optical gap of ~0.6 
eV. Based upon comparisons with chemically similar 
nicklates as well as existing electronic calculations^^ 
we assign the excitations centered at 1.5 eV as on-site 
nickel d to d transitions, likely activated by modest hy- 
bridization with the coordinating oxygens. These exci- 
tations take place in both majority and minority chan- 
nels according to recent DFT calculations^ The features 
above 3 eV are assigned as O p to Ni d charge transfer 
excitations. The energy range of our magneto-optical in- 
vestigation is also shown in Fig. \7\ providing a preview 
of the physical origin of the field-induced spectroscopic 
changes, discussed below. 

Figure Efa) displays the magneto-optical response of 
{NI4M012} as a function of magnetic field from to 32 T 
at 4.2 K. These experiments were carried out on the same 
pressed powder sample using a grating spectrometer (0.8 
- 3.5 eV) equipped with InGaAs and CCD detectors and a 
33 T steady field resistive magnet at the NHMFL in Tal- 
lahassee, FLi2i The reflectance ratio, R(B)/R(B = 0), is 
a normalized response and highlights changes in the opti- 
cal properties with applied magnetic field. With increas- 
ing field, the reflectance of {N14M012} decreases by ~2%. 
It is notable that this effect occurs in the visible spectral 
range, hence the name "magnetochromism" . Based upon 
the aforementioned peak assignments in the optical con- 
ductivity spectrum in Fig. [7] and the a\{B)/(j\{B = 
T) ratio data in the inset of Fig. 02b), we attribute the 
observed magnetochromic effect to a field-induced modi- 
fication of the Ni d to d on-site excitation. Distortion of 
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FIG. 8: (color online) (a) The normalized magneto-optical 
response, R(B)/R(B = T), of {M4M012 } pressed powder in 
an applied magnetic field from to 32 T at 4.2 K. 1 T steps 
are shown. No hysteresis is observed, (b) Absolute value of 
the integrated area of the magneto-optical reflectance ratio 
feature as a function of applied magnetic field (solid symbol) . 
Results from the integrated optical conductivity data are sim- 
ilar. The inset shows the change in the optical conductivity of 
{N14M012} with magnetic field, <ti(B = 30 T)/cri(B = T). 
The B = T curve ratios data taken before and after the field 
sweep, giving an indication of overall spectral reproducibility. 



the pseudo-octahedral Ni 11 crystal field environment is a 
plausible driver. Note that this is a local, molecular-level 
distortion rather than a bulk effect. No field- induced 
changes were observed on the leading edge of the O p to 
Ni d charge transfer bands above 3 eV. 

We quantify the magneto-optical effect in {N14M012} 
by plotting the absolute value of the integrated intensity 
of the reflectance ratio as a function of applied magnetic 
field (Fig.|^b)); — The magneto-optical effect is small at 
low fields, becomes appreciable above 10 T, and continues 
to grow for B > 30 T. There is no evidence of saturation 
to 33 T. The overall rising magnetochromic response can 
be fit with several different functions including a cubic 
polynomial in \B\ and a simple exponential, suggesting 
a likely functional form for the field-dependence of the 
magnetic parameters in {M4M012}, described below. We 
propose that the applied magnetic field interacts with 
the spin centers, deforming the local structure around 
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the Ni 11 sites. This process modifies the crystal field en- 
vironment, the result of which is a field-induced change 
in electronic structure (Fig. |H1 (a)). Thus, these spec- 
troscopic results support the picture of field-dependent 
exchange and anisotropy terms in the spin Hamiltonian 
of {M4M012} that derive from magnetoelastic (and con- 
sequent spin-orbit) interactions. Magneto-optical effects 
due to field-induced changes in local structure have also 
been observed in other materials including (CPA)2CuBr2 
and K 2 V 3 8 ^^ 



uniform spacing of the crossing fields has not been forth- 
coming. Assuming a Hamiltonian with two exchange cou- 
plings and anisotropic as well as biquadratic terms as 
given in did not result in a satisfactory description of 
all magnetic observables on a common footing. 27 There- 
fore, we extended our model to allow J and D to vary 
with applied magnetic field. The possible field-induced 
distortion of the crystal field environment around the Ni 
centers motivates this ansatz. 



IV. THEORETICAL MODELS 



FIELD-DEPENDENT MODEL 
PARAMETERS 



In order to understand the magnetic properties of 
{M4M012} we adopt the following general Hamiltonian 
that has been used for other Ni n -compounds, 16,18 



H = Hn + H ani + H biq + H z , where 



-Hani — D 



U 

2 juv (s(u) • s(v) 



H z = gii B B-Y H u ) 



(3) 
(4) 

(5) 

(6) 
(7) 



Here Hn denotes the Heisenberg Hamiltonian and H ani 

describes the single-site anisotropic (ligand field) con- 
tribution, which is compatible with the approximate 
tetrahedral symmetry of {Ni4Moi 2 }. The unit vector 
e r (u), which points radially outwards, serves as a local 
anisotropy axis for site u. The term 8/3 is convenient 
in order to render the Hamiltonian traceless. -ffbiq rep- 
resents biquadratic terms, which are the next higher or- 
der compared to the Heisenberg Hamiltonian^ 6 . and Hz 
is the Zeeman term. We employ a single spectroscopic 
splitting factor g since the g-tensor anisotropy was found 
to be very small for the present system; a similar result 
has been found for Ni 11 squares^ We also assume that 
a possible anisotropic exchange between the Ni 11 centers 
is small because the orbital contribution to the ground 
state is small^ 

In the following sections, we simplify J uv = J 
and j uv = j, or, if we use two different con- 
stants, J' = J12 = J13 = J14, J = J23 = J24 = J34 and 
f = 3l2 = J13 = 3x4 > 3 = 323 = 324 = 33A, as illustrated 
in Fig. |2| We also account for impurity effects (addi- 
tional paramagnetic Ni 11 ions) and their batch-to-batch 
variation by adding a paramagnetic term to the Hamil- 
tonian ©. 

We have made numerous attempts to model the exper- 
imental magnetization curve of {Ni4Moi2} (Fig. with 
field-independent values of J and D, as described previ- 
ously. Despite these efforts, an explanation for the non- 



A. Low-field properties of {Ni4Moi2} 

It is not a priori clear how the model parameters 
should depend on the magnetic field since this depen- 
dence is indirect. The exchange couplings, the bi- 
quadratic contributions as well as the anisotropy result 
from the electronic structure, and are (complicated) func- 
tions of orbital overlaps, lattice stiffness, and spin-orbit 
coupling. If the electronic structure is noticeably altered 
by an applied field it is probable that the parameters en- 
tering Hamiltonian ||2Jl are also changed, but very likely 
in a highly non- linear manner. To the best of our knowl- 
edge, so far such dependencies were only observed in cer- 
tain Mn grid molecules^ where the anisotropy changes 
sign at a level crossing. It appears plausible to us that 
besides continuous variations, the parameters can also 
change abruptly since the molecular structure may relax 
into new ground states at certain field values. Therefore, 
one can anticipate that the field-dependence of the model 
parameters might only be piecewise analytic between the 
abrupt changes. In this section we therefore investigate 
simple piecewise parameterizations of the Hamiltonian. 
We also neglect a possible anisotropic field dependence 
of the parameters, i.e. the exchange parameters of dif- 
ferent bonds would be modified differently, and therefore 
the effect should depend on the relative orientation of the 
molecule with respect to the field. 

It turns out that the low-field susceptibility versus tem- 
perature data, shown in Fig. can be modeled by two 
sets of parameters. The magnetization of two differ- 
ent samples of {Ni4Moi2} at B = 0.5 T can be un- 
derstood assuming J/k B = —3.2 K, J'/ks = —3.2 K, 
j/k B = 1.6 K, j' = 0, D/k B = -1.0 K together with 
the separate term to account for the paramagnetic impu- 
rities (free Ni 11 ions) contained in the samples (parame- 
terization 1). The parameters of the Hamiltonian com- 
pare nicely to those of 2 x 2-Ni 11 grid molecules; 18 The 
sign of the anisotropy as well as that of the biquadratic 
term signal the same behavior, only the absolute value of 
the biquadratic term is somewhat larger. The resulting 
zero-field splitting of the triplet is only 0.15 K, in very 
good agreement with our EPR results. The fact that 
the biquadratic term is one order of magnitude larger 
than usually observed (j oc J/lOO)! 8 ^*^ signals that 
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T (K) 

FIG. 9: (color online) Low-field susceptibility: The mea- 
sured data are given by symbols. The fits for B — 0.5 T 
are obtained using J/k B = —3.2 K, J' /fee = —3.1 K, 
j/k B = 1.6 K, j'/k B = 0, D/k B = -1.0 K. The data for 
field values of B = 3.5 T up to 5 T are approximated using 
J/K = -3.2 K, J'/k B = -3.2 K, j/k B = 1.5 K, j'/k B = 0, 
and D/k B — —3.2 K. g — 2.195 in all cases. The estimated 
impurity concentrations are 0.216 and 0.146 individual Ni 11 
ions per {M4M012} molecule for sample 1 and 2, respectively. 
The simulated data are averaged over 100 orientations. 



change. Such a change could be caused by variations 
of the bond distances and angles; the latter are known 
to have dramatic effects on the exchange parameters. 
Based on the mag neto-structural correlations for Ni-0 
cubanes in Ref. 28, a change of the exchange parameter in 
{N14M012} by the necessary amount would correspond to 
a change of the bond angles by approximately one quar- 
ter of a degree if the effect was attributed solely to bond 
angle modifications. A literature survey of chemically- 
similar Ni(II) cluster compounda 16 i 18 i 20 i 21 i 22 i 28 i 29 i 30 re- 
veals that J and D vary from at least -10 to 17.5 cm -1 
(-7 to 12.2 K) and -5.5 to 0.5 cm" 1 (-3.8 to 8.5 K), re- 
spectively. In these systems, bond angle is the govern- 
ing structural parameter due to the intimate relationship 
between angle and magnetic orbital overlap^ There is 
a linear correlation between the value of J and Ni-O- 
Ni and Ni-O-O-Ni bond angles in several nickel clus- 
ters. Although the transition metal centers in these nickel 
clusters have different local environments, the aforemen- 
tioned range of J and D represents the variation in phys- 
ical parameters that can be accessed by chemical tuning 
in these compounds. Magnetic field-induced effects are 
anticipated to be in this range as well. 



the bonds in the basal triangle of the molecule might 
be rather soft. The other four susceptibility measure- 
ments, Fig. 03 which are obtained from B = 3.5 T to 5 T, 
can be better approximated by a second set of parame- 
ters: J/k B = -3.2 K, J'/k B = -3.1 K, j/k B = 1.5 K, 
j'/kB = 0, and D/k B = —3.2 K (parameterization 2). 
These parameters differ from those for B = 0.5 T in two 
ways: (1) the biquadratic contribution is 7 % smaller, (2) 
the anisotropy coefficient is three times larger. 



B. High-field properties of {NI4M012} 
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Although the above parameterization describes the 
low-field data up to B w 5 T reasonably well, it fails to 
reproduce the magnetization data at higher fields. The 
dramatic increase of the field spacings between adjacent 
magnetization steps at 4.5, 8.9, 20.1, and 32 T cannot 
be explained by small changes of the anisotropy D or 
the coupling J. The third field spacing is about twice 
the first one, which using relation (0 suggests an ex- 
change coupling that is about 1.25 of the original one. In 
the following we therefore investigate a model where only 
the exchange parameters are allowed to depend on field. 
This is of course a simplification since all parameters of 
the Hamiltonian (J, j, D, and g) should be modified in a 
varying field. For our purpose we assume a hypothetical 
exponential dependence of J on the absolute value of the 
external magnetic field 

J(B) = J exp(^j , ( 8 ) 

which we motivate by a possible change of the over- 
lap of those orbitals that are involved in the superex- 



FIG. 10: (color online) Magnetization of sample 2: Experi- 
mental data (NHMFL) are given by squares. The solid curve 
provides the best fit using a Hamiltonian with an exponen- 
tially field- dependent coupling. The simulated data are aver- 
aged over 100 orientations. The inset shows the experimental 
differential magnetization AM/AB as data points as well as 
the theoretical AM /AB (solid curve) using the same parame- 
ters as above. 



The best fit using two exchange parameters was ob- 
tained for J/k B = -4.2 K * exp(B/96 T), J'/k B = 
-3.2 K * exp(£/52 T), j/k B = 0.16 K, j'/k B = 0.39 K, 
D/k B = —8.9 K, and g = 2.195 together with the para- 
magnetic impurity contribution already determined for 
sample 2 (parameterization 3). This fit is shown in 
Fig. EI] Except for small fields, where the two afore- 
mentioned parameter sets are appropriate, this model 
provides a rather good description of the magnetization 
over a large field range. The assumption of an exponen- 
tial dependence is not essential; a polynomial dependence 
yields similar results. 
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EPR revisited 



The level scheme and possible EPR transitions can be 
shown more closely using the parameters of the general 
Hamiltonian 




FIG. 11: (color online) The figure shows the Zeeman level 
splittings for the realistic Heisenberg model © (parametriza- 
tion 2) for several orientations of the field relative to the 
molecules. The various curves, which differ especially at 
avoided level crossings, represent relative angles of 10° . 



three transitions, which are shown in Fig. ^]for a fre- 
quency of v = 190 GHz and a temperature of T = 4.2 K. 
The transition PI, which occurs between energy levels 
marked by filled circles in Fig. 1111 is rather narrow; this 
is connected to the fact that variations of the field di- 
rection do not lead to strong variations of the respective 
levels. This also implies that the mixing of S z eigenstates 
is weak and thus the amplitude of this forbidden transi- 
tion is small. P3 on the contrary is an allowed transition 
which dominates the spectrum. In the investigated field 
range it is given mainly by the transition between levels 
marked by open squares in Fig. ^2 but transitions be- 
tween other (unmarked) levels also contribute, he rather 
broad line can in part be explained by the rather strong 
variation of the upper marked level with variations of the 
field direction. The transition P2 between levels marked 
by open triangles in Fig. ^] is stronger than PI but also 
broader, which can again be understood by looking at 
the variation of levels. 



VI. CONCLUSIONS 



In Fig. ^2 the Zeeman split levels are presented for 
fields up to 8 T using parametrization 2. The thick curves 
show levels for an orientation of the field axis pointing 
from the mid point of the basal triangle through the top 
Ni center (Ni ion 1 in Fig. |2J) . The thinner curves show 
the levels for orientations with relative angle steps of 10° 
along a great circle through one of the basal Ni centers. 
The variation of the level positions with orientation for 
a fixed absolute value of the external field is especially 
large at avoided level crossings. Although the final EPR 
line is given by averaging over all orientations one can 
already deduce from Fig. ^2 that the EPR measurement 
is in full agreement with the microscopic Hamiltonian. 
The attached symbols for the three observed transitions 
match the calculated level spacings very well. 
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FIG. 12: Relative EPR transmission observed at v = 
190 GHz, which corresponds to an energy level difference 
of 9.12 K.: The dominant (allowed) transition P3 is rather 
broad, whereas the the weaker (forbidden) transition PI re- 
sults in a much sharper peak. The other weak (forbidden) 
transition P2, which is masked by P3, is broader than PI. 



To some extent one can even explain the widths of the 



We have presented a comprehensive experimental 
and theoretical investigation of the magnetic molecule 
{M4M012}. We find that the main model parameters 
of the Hamiltonian (i.e. exchange and anisotropy pa- 
rameters) have a strong dependence on magnetic field, 
an effect that may be accompanied by molecular mag- 
netostriction. All of our efforts to avoid this conclusion 
lead to a blatant contradiction between theory and ex- 
periment. The discovery of a dependence of Hamiltonian 
parameters on field in {M4M012} is quite novel. Nev- 
ertheless, it may be a general characteristic of magnetic 
cluster-based materials with strong lattice coupling. It 
is well known that most materials (including magnetic 
molecules) are not j-igi^^^JS^L^ In addition, only rel- 
atively few high-field magnetization studies have been 
performed to date. Thus, it is of interest to discover 
whether other magnetic molecules might display a simi- 
larly strong dependence of the model parameters on mag- 
netic field in the regime above 20 T. 
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